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Abstract 

The stationary state of stochastic processes such as reaction-diffusion systems can be related 
to the ground state of a suitably defined quantum Hamiltonian. Using this analogy we investi- 
gate the applicability of a real space renormalisation group approach, originally developped for 
quantum spin systems, to interacting particle systems. We apply the technique to an exactly 
solvable reaction-diffusion system and to the contact process (both in d = 1). In the former 
case, several exact results are recovered. For the contact process, surprisingly good estimates of 
critical parameters are obtained from a small-cell renormalisation. 
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1 Introduction 



Reaction-diffusion systems and other interacting particle systems are relevant for the description 
of several phenomena in physics, chemistry and biology Q. In the past, they have been mainly 
modelled by (non- linear) partial differential equations Q , a description which implicitly contains a 
mean-field assumption. Such a description is however no longer appropriate in low dimensions where 
fluctuations are important. To take these into account one turns to a description of the reaction- 
diffusion system in terms of a stochastic process. This can for example be realised by adding a 
noise term to the partial differential equation. In recent years however, particular attention has 
been paid to models defined on a lattice. It has been found that these can be related to a number of 
interesting topics in modern statistical mechanics such as growing interfaces H, phase transitions 
into an absorbing state 0], exactly solvable quantum spin chains ||, persistence exponents H and 
so on. 



The systems which we will study here are defined on a lattice but evolve in continuous time. At each 
site of the lattice, one can have (hard core) 'particles' which can perform a random walk and/or can 
undergo one or several 'reactions'. In this paper we will limit ourselves to systems with one type 
of particle . Each lattice site can then be either empty (0) or be occupied by a particle (A). As an 
example, consider a system in which particles perform random walks and where two particles on 
neighbouring sites can 'annihilate' (i.e. undergo the reaction A + A ^ + 0). In a simple mean field 
approach the density of particles c{t) in this system decays asymptotically as 1/t. It is common 
to introduce a critical exponent 9 which describes the decay of the density (c(t) ~ t~ e ) and which 
therefore in mean field theory equals 1. An exact solution of the diffusion-annihilation model in 
d = 1 (where it is equivalent to the T = Glauber dynamics of an Ising model on the dual lattice) 
shows however that = 1/2 0]. This latter value is also found experimentally in systems which are 
thought to be described by the diffusion-annihilation model. Moreover, 9 shows a large universality 
across materials and initial conditions. The same value for 9 is also found when the experimental 
situation is described more correctly by a diffusion-coagulation (A + A — > + A, A + A ^ A + 0) 
model f||. Hence, as in the theory of equilibrium critical phenomena one needs a scheme which 
at the same time explains the observed universality and gives precise values for critical exponents 
such as 9. One is therefore naturally led to a search for renormalisation group (RG) approaches to 
stochastic systems. 



Mathematically, the stochastic process is a continuous time Markov process. Using the formal 
equivalence between the master equation and the Schrodinger equation in imaginary time, one 
can set up a field theoretic formulation for stochastic systems which in turn can be used in the 
construction of a renormalisation approach Q. Critical exponents can then be calculated in an e- 
expansion around the upper critical dimension, which in this case equals 2. This approach is by now 
well established and has been applied to many interesting systems. As an example, we mention the 
much studied branching and annihilating random walks (BARW) fi~0|| , where besides diffusion and 
annihilation, particles can undergo branching processes (A — > (m + l)A). The competition between 
annihilation and branching leads to a non-equilibrium phase transition where the stationary state 
particle density c s t(= limt^ 00 c(t)) goes to zero in a continuous way as a function of the rates of 
the different processes. Near this transition, several critical exponents (static and dynamic) can 
be introduced. It has been shown, both numerically and using the field theoretic RG, that the 
universality class of the transition for BARW is completely determined by the parity of m [11]. For 
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m odd, the transition falls in the universality class of directed percolation, whereas for m even a new 
universality class appears. On the basis of very precise simulations, Jensen [jf^] conjectured that 
8 = 2/7 and (3 = 13/14 when m is even {(3 describes the way c s t goes to zero near the transition). 
So far, a precise analytical calculation of these exponents has however not been possible. Menyhard 
and Odor propose the value /3 = 1 on the basis of a perturbation around mean field theory [13|. 
Using a loop expansion at fixed dimension, Cardy and Tauber find f3 = 4/7 [|llj. Thus, there is 
clearly room for the introduction of new, analytical approaches to the BARW and related models. 



In this paper, we investigate the possibilities of a real space RG approach to interacting particle 
systems. Our starting point is again the equivalence between a stochastic system and a quantum 
mechanical one. In the past, several real space RG approaches to quantum lattice systems have 
been introduced. We must mention here as an example the density matrix renormalisation group 
(DMRG) which has been very succesful [14]. Currently, several research teams in the world are 
investigating the applicability of the DMRG to reaction-diffusion systems [fl5f , [plj] . Unfortunately, 
the name DMRG is a bit a misnomer, since one rarely calculates RG flows in this approach and hence 
it is not easy to decide on questions of universality using this technique. Moreover, the approach 
is purely numerical. Instead in this work, we will use an appropriate perturbative extension of 
the so called standard (also called SLAC) approach [17]. We have applied this RG method to the 
study of the stationary state properties of some interacting particle systems. We found that the 
method works surprisingly well in several cases which we studied. These include a solvable reaction 
diffusion model (with diffusion, coagulation and decoagulation) and the contact process, which is 
also thought to be in the directed percolation universality class. 



This paper is organised as follows. In section II we briefly introduce the description of stochastic 
processes in a quantum mechanical language. In section III, we discuss a real space approach to the 
ground state of quantum (spin or fermion) chains. Our method is perturbative and in first order 
coincides with the standard (or SLAC) approach. We discuss how the technique can be derived 
from more conventional (Niemeijer - Van Leeuwen) RG approaches to quantum lattice systems 
by taking the zero temperature limit. The method as such is not new but, in our opinion, not 
very well known. We also discuss how critical exponents for stochastic systems can be obtained 
from such an approach. In section IV we study properties of the stationary state of some solvable 
reaction-diffusion systems. In section V, we give our results for the contact process. Finally, in 
section VI we present some conclusions and an outlook on further applications of the real space 
RG technique. 



2 Quantum formalism of reaction-diffusion systems 



In this section we discuss the relation between stochastic processes in continuous time on the one 

This relation has been (independently) 



hand and quantum mechanics on the other hand [18|,[|1J 
discovered by several authors, and is now well established. We only give a brief overview, with the 
main purpose of fixing our notation. 



Consider a one dimensional lattice of L sites and let r\ = {7/1, . . . , rji} = 0(1) when no (a) 
particle is present at site i) be the microscopic configuration of the particle system. Furthermore, 
we denote by P(r], t) the probability that the system is in configuration r/ at time t. The time 
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evolution of P(j],t) is determined by the transition rates wq, w±, ... ,w s of the model. In general 
one of these is used to fix the time scale (say wq = 1). We will collect the remaining rates in a 
vector w = (w\, . . . ,w s ). 



Next, we introduce a state vector \P(t)) 

\P(t))=J2P(v,t)\r)) (1) 
v 

The time evolution of \P(t)) is described by the master equation 

d\P(t)) _ 



dt 



H\P(t)) (2) 



where H is a (2 L x 2 L ) matrix which depends on all the transition rates of the system. In particular 
the matrix element — i?f £ > (for C equals the transition rate to go from configuration £ to 
configuration £. Due to the conservation of probability the diagonal elements are given by 

#C,C = " E H U ( 3 ) 

This relation implies that the sum of the elements in a column of H equals zero, a condition which 
we will refer to as the stochasticity condition. 



It is now common to consider (^j) as a Schrodinger equation (in imaginary time) and call H the 
Hamiltonian. In contrast to the quantum mechanical situation H is not necessarily hermitian. 
Hence, left and right eigenvectors are in generally not related by simple transposition. In situations 
where there are only local interactions we can write H in terms of local operators H = J2i Hi. For 
reaction-diffusion systems with one type of particle, it is common to interpret the model as a spin 
model (where a spin up (down) corresponds the absence (presence) of a particle) so that one writes 
H in terms of Pauli- matrices or in terms of the matrices E a/3 where (E a ^) aiT = 5 aa 5f3 T . 

The stochasticity condition implies the existence of at least one eigenvalue which equals zero and 
an associated left eigenvector (s\ given by 

(s\ = 5>| (4) 

v 

Moreover, it can be shown that the real part of all the eigenvalues is non-negative. 



The formal solution of the master equation @ is 

\P(t)) = e- Ht \P(0)} (5) 

The expectation value of an operator, say X, corresponding with a physical quantity (such as the 
density of particles) can then be written as 

(x(t)) = J2 x WP(v,t) 

v 

= (s\X\P(t)) 

= (s\Xe~ Ht \P(0)) (6) 
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In this paper, we will be mainly interested in stationary state properties of the stochastic system. 
In that limit, (||) is determined by the groundstate(s) |0) of H. In particular, (|J) becomes 

X st = lim (X(t)} = (s\X\0) (7) 

t— >oo 

in the case of a non-degenerate ground state |0). 



In summary, when we want to study properties in the stationary state of a stochastic system we 
are in need of a (real space) RG which is suitable for the ground state of quantum spin systems. 
In the next section, such an approach will be presented. 



3 Ground state renormalisation for quantum spin systems 



It is convenient to remind the reader briefly on the history of real space RG approaches to quantum 
systems. In a first development, approaches working at finite temperature were introduced pQ],pl]]. 
These were direct extensions of the first real space RG developped by Niemeijer and Van Leeuwen 



[22]. The SLAC approach was later introduced as an independent approach to the ground state of 
quantum lattice systems. It was at first claimed to be non-perturbative. In a little known paper 
p3| , it was however shown that the SLAC approach can be obtained from the finite temperature 
technique by applying a suitable perturbation expansion (see below) to first order and by sending 
the temperature T — > 0. We now briefly describe the two methods and their relation. As will 
be explained below, the possibility to turn the SLAC approach into a perturbative technique is of 
importance for the study of some stochastic systems. 



The SLAC or standard approach was introduced in |Tfl| . It was used a lot in the early eighties to 
study ground state properties of several quantum spin and fermion chains (for a review, see p4|]). 
As usual in real space RG methods, the lattice is divided into cells, each containing b sites. The 
Hamiltonian H is then divided into an intracell part Hq and an intercell part V. If a labels the 
cells, one can in the particular case of one dimension write 

H = Y J (H , a + V a>a+1 ) (8) 

a 

As a first step, the Hamiltonian within one cell Hq i(X is diagonalised exactly. Let i?o,a have eigen- 
values E n ^ a with corresponding right (left) eigenvectors |n) Q ( Q (n|). One then selects two low lying 
eigenstates (for example the ground state |0) Q and the first excited state |1) Q ) and considers them 
as effective spin states for the cell: \+ a )' = |0) Q and \— a )' = Renormalised lattice states \a)' 
can then be constructed by making tensor products over all cells: \a)' = ® a \cr a )'. These states 
span a 2 L l h dimensional vector space W. 



The renormalisation transformation, which always involves an elimination of degrees of freedom, is 
now performed by projecting the original Hamiltonian onto W. Mathematically this is achieved by 
introducing a 2 L x 2 Z/ / b -matrix T2 whose columns contain the vectors \a)' together with a 1 L l h x 2 L - 
matrix T\ whose rows contain the vectors ' (a\ (which are constructed from the left eigenvectors of 
Hq^). Then the renormalised Hamiltonian H' = TZ(H) is defined as 

H' = K(H) = TiHT 2 (9) 
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The transformation (||) will define a mapping in the parameter space w' = f(w), from which fixed 
points and critical exponents can be determined as we will explain further below. 



The procedure which was defined above at first sight seems to be rather ad hoc and non-perturbative. 
Further insight in the method was obtained when its relation with more conventional real space 



RG approaches was discovered [23]. These approaches work at finite temperatures and are a direct 
extension of the Niemeijer-Van Leeuwen real space approach to the case of quantum spin systems. 
In such an approach, the eigenstates of i?o,o are divided into two groups according to some criterion 
(for example, a majority rule). Each group is associated with an effective spin state for the cell. 
Within each group there are 2 b ~ 1 states. Hence, cell states can be denoted as \a a , q a ), where 
a a = ±1, q a = 1, . . . , 2 . Tensor products of these states will be denoted by \a, q) = ® a \(T a , q a ). 
Finally let 7i = —H/(kT) be the reduced Hamiltonian. The matrix elements of the renormalised 
(reduced) Hamiltonian TC are then obtained by performing a partial trace 

'{a\e H '\a'Y = Tr q (a, q\e n \a' , q) (10) 



Usually, it is impossible to work out this mapping exactly. In any explicit calculation, it is therefore 
necessary to perform an expansion in the intercell Hamiltonian V . This is achieved by using the 
Feynman identity 



- X7io Ve m °d\ 



(11) 



where T\ is a time ordering operator and 7i$ and V are respectively the reduced intracell and 



intercell Hamiltonian. A similar expansion is made on the left hand side of (|i0[). In reference [23] it 
was shown that if one expands ( JTT| ) to first order in V, and then takes the limit T — > 0, one recovers 
the SLAC approach. This relation is useful for several reasons. First and most importantly, it shows 
that the SLAC approach is perturbative and one obtains a procedure on how to calculate higher 
order corrections. Secondly, since the finite temperature RG is constructed to conserve the partition 
function, one is garanteed that the SLAC approach conserves the ground state energy. Finally, this 
relation allows a consistent study of finite and zero temperature properties of a quantum system 
in thermal equilibrium. 



Going back to reaction-diffusion systems, which are non-equilibrium systems, we are thus garanteed 
that our RG approach conserves the ground state energy. In principle, we also have a recipe to 
calculate perturbative corrections to the SLAC approach. In this paper, we will limit ourselves to 
calculations in first order. Higher order corrections usually lead to a proliferation of terms in the 
Hamiltonian. However, as we will discuss in our conclusions, the existence of these higher order 
terms will be necessary for a proper study of the BARW with m = 2, or for any other model in 
which H contains terms involving three or more sites. 

While in this paper we will only study stationary state properties (corresponding to ground state 
properties of H ) , it is our feeling that the finite temperature extension of the SLAC-approach will 
be useful to study finite t properties of the reaction-diffusion systems. 



We now turn to a discussion of how stationary state properties can be determined from the RG 
mapping w' = f{w). To fix ideas, let us assume that this equation has a non-trivial fixed point at 
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w*, with one relevant scaling field (which in linear approximation is proportional to Awi 



whose scaling dimension is y Wl . From standard RG theory it then follows that near 



w - 

criticality, the correlation length £ will diverge as £ ~ lAwil^- 1 - with i/± = l/y Wl . 

In general it will be so that after renormalisation w' ^ 1. Hence time needs to be rescaled again, 
which is achieved by dividing H' by w' . The number w'o(w*) therefore teaches us how time rescales 
under a rescaling of space. We can use this quantity to calculate the exponent z as 



(12) 



From v± and z, scaling [25] gives us the exponent u» = v±/z which determines the divergence of 
the relaxation time near the critical point. Next, we turn to the calculation of the particle density 
c s t- For a translationally invariant system, we can, using (Q), write c s t as 



which under the RG transforms 



c st = <s|£ n |0) 
transforms as 



c st = (slT^T^TiT^O) 

= '(s\(E n y\oy 



(13) 
(14) 



The simplest possible case is the one in which the renormalisation of the operator E 11 does not 
involve other operators so that (E 11 )' = T\E ll T2 = a(w )E n . In that case ( |l4]) becomes (where we 
now explicitly denote the dependence of c s t on the transition rates) 



c s t{w) = a(w)c st (w') 



(15) 



This relation can be iterated along the RG-flow, and hence the density of particles can be obtained 
as an infinite product if one knows the density at the fixed point Wt which attracts w (where the t 
reflects the fact that this attractive fixed point is trivial and not critical). 



c s t{w) 



n 

Li=0 



(16) 



In principle correlation functions can be calculated in a similar way. 



Near the critical fixed point w*, we get from (15) for the singular part of c s t to leading order in 



c st (A Wl ) = a{w*)c st (b y ^ Awi) 



(17) 



From this relation we see that a(w*) determines the rescaling of the particle density at criticality. 
We write a(w*) = b D ~ d where D can then be interpreted as the fractal dimension of the sites 
that are occupied by particles at criticality. Finally, from (|l7]) we get the behaviour of the particle 
density close to the critical point as 

c st (Awi) ~ (Awi) {d ~ Dyy ^ (18) 

so that we obtain the scaling relation (3 = (d — D)v±. Finally, 9 can be obtained from the scaling 



relation 6 = f3/i>\\ p| 
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It is straightforward to extend these reasonings to the case where the transformation of E is more 
complicated. 

We thus see that a complete characterisation of the stationary state particle density and of all the 
critical properties of the system can be obtained from our RG approach. In the next two sections we 
test our method on simple stochastic systems. The first one is an exactly solvable reaction-diffusion 
system with a trivial transition. The second one is the well known contact process which has a non 
trivial transition thought to be in the DP universality class. 



4 Renormalisation for a simple reaction-diffusion system 



We consider a model with diffusion, decoagulation (which is the process A + % — > A + A, ® + A — > 
A + A) and coagulation. We will use the diffusion rate to fix the time scale, so that our model has 
two independent rates which we will denote as w\ (decoagulation) and W2 (coagulation). The local 
Hamiltonian Hi in this case can most conveniently be written as a 4 x 4-matrix 



H 










-1 

— W\ 




-1 

1 + W\ 

— Wl 



2w 2 J 



or in terms of the matrices E a @ as 



(19) 



- ErE^-E^E^-MErElU+E^E^-w^E^ + E^E^ (20) 
With a similarity transformation this Hamiltonian can be mapped onto that of a free fermion 



system [26], from which many exact results can be obtained. 



To renormalise this model, we take b = 2, so that i?o,o is exactly given by (pop. The spectrum of Hi 
can be calculated trivially. The ground state is doubly degenerate and we identify the corresponding 
eigenstates as effective spin states. In particular, we have (we use spin language where |+) (|— )) 
denotes a vacancy (an A-particle)) 



1+)' 
'(+1 

I-)' 

'(-I 



-) + [-+) +r| 

+ (- + ! + (-- 



(21) 



where r = wi/w 2 and N = 2 + r. Notice that we normalise the states as probability vectors (and 
not quantum mechanically). 



Next, we calculate H' using (^). We find that the renormalised Hamiltonian contains the same terms 
as (^0|). Moreover, the transformation conserves the stochasticity condition. The renormalised 
diffusion rate is 1/N. We divide the Hamiltonian by this factor (rescaling of time). Then H' has 
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completely the same form as H but with renormalised values for w\ and w 2 • The RG equations for 
w\ and w 2 are 

, 1 + wi + w 2 . . 

W] = w i (22) 

^2 

, 1 + ^1+^2 , OQ , 

w 2 = w 2 — ■ (23) 

2w 2 + w\ 

The flow generated by these equations is shown in figure 1. There is a fixed point at w 2 = 1, w\ = 
(pure coagulation fixed point). The line w 2 = 1 is an invariant line. The RG equations, linearised 
at the fixed point, have the relevant eigenvalue y Wl = 1, together with a marginal eigenvalue in the 
w 2 direction. 



Next, we turn to a calculation of particle density in the stationary state using the scheme outlined 
in the previous section. Projecting E 11 onto W we find the simple recursion 

(E 11 )' = (24) 

From this we have 

1 + w\jw 2 



a(wi,w 2 ) 



2 + wi/w 2 



Since also w' x jw' 2 depends only on the ratio w\/w 2 we arrive at the conclusion that the particle 
density only depends on that ratio. Its precise value can then be calculated by making a product 
of a(wi,w 2 ) along the RG flow. From figure 1 we see that all points with w\ > flow to points 
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Figure 2: Particle density in the stationary state of the diffusion, coagulation, decoagulation model. 
The full line represents the exact result, the circles are the RG results. 

where w\/w2 — * oo. In that limit c s t = 1. Then, using (|i~6|) we can obtain the particle density. In 



figure 2 we plot the result for function of r together with the exact result |26j 



c s t(r) 



r + 1 



(25) 



Within the numerical accuracy, both results coincide. Hence it seems that we recover the exact 
result! This might seem surprising since our calculation is only precise to first order in V. With 
hindsight the accuracy of our result can be understood by the fact that the ground state of the 
whole system is a product of one particle states. Nevertheless, this calculation illustrates nicely the 
use of the method. Moreover it gives an RG-flow for the problem which cannot be obtained in any 
other way. 



To conclude this section, we calculate the exponent (3 which determines the behaviour of c s t near 
r = 0. From the exact result (p5|), one obviously has (3 = 1. On the other hand, from (pip, we have 
a(0, 1) = 1/2, hence D = and using y Wl = 1 we recover the exact result (3 = 1. 



At this place it is appropriate to mention that for other simple reaction diffusion systems we can 
also recover exactly known results. An example is a model with diffusion, decoagulation and death 
(which is the process A + — > + 0, + ^4 — ^0 + 0). This model undergoes a first order transition 
|p6| . When the decoagulation rate is greater then the death rate, c s t = 1, whereas when the 
opposite inequality holds, one has c s t = 0. The RG recovers this exact result. As usual the first 
order transition is controlled by a discontinuity fixed point. 



Also for a model with diffusion, annihilation and pair creation + — > A + A, we recover the exact 
result, first obtained by Z. Racz pTfl, that c s t grows with the square root of the pair creation rate. 
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5 The contact proces 



The contact proces was originally introduced as a simple model for an epidemic ]2q| . In that 
interpretation a particle corresponds to sick person, and a vacancy to a healthy individual. In the 
process particles can disappear (A — > 0) with a rate wq = 1. Empty sites can become occupied 
with a rate Xz/2, where z is the number of occupied neighbours of the empty site (this represents 
contamination in epidemic terms). This model cannot be solved exactly, but its critical exponents 
are known to high accuracy ||. On the basis of numerical data, and from symmetry arguments, it 
is generally believed that the model is in the DP-universality class p9|] , [|30|] . 

The quantum Hamiltonian corresponding to this model contains one-site terms (for the process 
A — > 0) and two-site terms (for the contamination process). The local Hamiltonian for the model 
is 



Hi 








-1 


°\ 





A/2 





-1 








l + A/2 







-A/2 


-A/2 


1 / 



(26) 



We have split the whole Hamiltonian in such a way that Hi contains one two-site contribution and 
one one-site term. In this way, when performing the RG , both Hq )01 en V aa +i will contain an equal 
number of one-site and two-site terms. This garantees, at least for the contact process, that the 
ground state of the intracell Hamiltonian is a doublet. These states are then the natural candidates 
to be used as effective cell spins. 

We will perform a renormalisation for this model using a cell with 6 = 3. The calculation is 
straightforward and can most effectively be done using ^Mathematica. It came as a surprise to 
us that also in this case, there is no proliferation of interactions in the renormalised Hamiltonian. 
The RG-equation for A is 

= A 3 (2 + A)(8 + 10A + 4A 2 + A 3 ) 

4(16 + 40A + 37A 2 + 18A 3 + 4A 4 ) 1 ' 

This equation has one non-trivial repulsive fixed point at A = A* = 3.22319. This value is sur- 
prisingly close to the best known numerical value for the contact process which is A c = 3.2978 (all 
numerical results are taken from fl3lj| ). From linearising the RG-equations near the fixed point, 
we obtain y Wl = .8886, from which we obtain v± = 1.1253, to be compared with the numerically 
determined value of 1.0972. As explained in section 3, one can obtain the exponent z from the 
rescaling of the unit of time. In this case we have 

, _ 4(16 + 40A + 37A 2 + 18A 3 + 4A 4 ) 
W ° ~ (8 + 10A + 4A 2 + A3)2 ^ 

from which we obtain z = .6858 (to be compared with z = .636). Finally we need to determine D. 
As explained in section 3, one therefore has to renormalise the operator E 11 . In a cell with b > 3 
this can be done in several ways. Either one projects the operator where m is a site at or near 
the middle of the cell. Alternatively we can take an 'average' operator E^ 1 
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Figure 3: Particle density in the stationary state of the contact process as obtained from the RG 
approach. 

We have performed both calculations. They give respectively D = .6391 and D = .6940, from 
which we obtain (5 = .4061, respectively [3 = .3444. These values should be compared with the 
precise numerical estimate (5 = .2769. Other exponents can be obtained using scaling relations. We 
get u\\ = 1.6409,0 = .2099 whereas the best known values are uu = 1.736,0 = .1597. We thus see 
that, taking into account the smallness of the cell considered, our estimates of A c and all critical 
exponents are close to the known values. 

In figure 3 we finally plot c s t(X) as obtained from our RG approach. At the scale of the figure it 
almost coincides with the results from other approaches. 

We are currently extending our calculations for the contact process to larger cell sizes. We hope 
that our approach, combined with suitable extrapolation techniques, is able to give very precise 
estimates of critical exponents. 



6 Conclusions 

In this paper we have investigated the applicability to reaction-diffusion processes of a quantum 
real space RG method. We first studied simple processes for which exact information is available. 
In three cases, these exact results are reproduced in calculations done on very small cells. 

So far, we encountered one system in which the RG predicted wrong results. This is a model with 
diffusion, coagulation and birth (0 + 0^0 + A, + — > ^4 + 0). When the diffusion rate equals the 
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coagulation rate, this model belongs to a class which is, at least partially, integrable [32]. For small 

1/3 

birth rates 103, c s t is known to grow as w 3 . The application of our technique to that model raised 
several difficulties. The ground state of -Ho, a is a singlet, and for the first excited state there is a 
crossing of energy levels (in finite systems). Hence, there is no obvious choice for the effective cell 
states. For a specific choice we made, the renormalised Hamiltonian H' turns out to have a form 
different from the original Hamiltonian H, but there is a similarity mapping this H' onto H with 
renormalised couplings, and an extra decoagulation term. If we define the full RG as the projection 
(|)|), followed by the similarity we are able to obtain a flow in parameter space. Unfortunately, our 

1 /2 

results indicate that for cells with 6 = 2 and b = 3, c s t grows as w 3 . We hope to clarify the RG 
for this model in the future. 



For the contact process, which contains a non-trivial critical point, rather accurate estimates for 
the location of the critical point and for critical exponents are obtained from a calculation on a 
cell of 3 sites. This justifies the hope that by going to larger cells and using good extrapolation 
techniques very precise exponent estimates can be obtained. We are currently performing such 
calculations. 



Another project we hope to carry through is a study of the BARW with m = 2. As stated in the 
introduction, there are few reliable analytical results on the critical properties of this system. A 
particularly nice model that is known to be in this universality class was introduced by Menyhard. 
It is a non-equilibrium Ising model (NEKIM) with Glauber dynamics at zero temperature and 
Kawasaki dynamics at infinite temperature p3|. It was recently shown that this model is selfdual 



1 34]. The ground state energy of the quantum Hamiltonian corresponding with this model is for all 
finite systems again doubly degenerated, so that effective cell states can be defined unambiguously. 
The quantum Hamiltonian for the NEKIM contains a three-site interaction term. Under the RG, 
performed to first order, such a term will be mapped onto a two-cell interaction term. However, by 
extending the SLAC-approach to second order in V, as discussed in section 3, one could generate a 
renormalised three-cell interaction. In fact, from a mathematical point of view, the Hamiltonian of 
the NEKIM is rather similar to that of a transverse Ising model with three spin interactions , Q] . 
That model also has a selfduality and was renormalised successfully by using our approach to 
second order |37| . An alternative model in the universality class of the even-m-B ARWis a recently 
introduced variant of the contact process in which particles disappear and are contaminated in 
pairs [38|. Since this model doesn't involve any diffusion, it may be more simple to analyse. 



It has of course to be admitted that real space RG methods in general involve some ill understood 
approximations. Nevertheless, it is our opinion that the results presented here give considerable 
hope that our RG can succesfully be used to further understand the critical behaviour of non 
exactly solved systems such as the BARW. 
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